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Abstract 

A half-tree is an edge configuration whose superimposition with a perfect matching 
is a tree. In this paper, we prove a half-tree theorem for the Pfaffian principal minors of 
a skew-symmetric matrix whose column sum is zero; introducing an explicit algorithm, 
we fully characterize half-trees involved. This question naturally arose in the context 
of statistical mechanics where we aimed at relating perfect matchings and trees on the 
same graph. As a consequence of the Pfaffian half-tree theorem, we obtain a refined 
version of the matrix-tree theorem in the case of skew-symmetric matrices, as well as 
a line-bundle version of this result. 



1 Introduction 

We prove a half-tree theorem for the Pfaffian principal minors of a skew-symmetric matrix 
whose column sum is zero. This is a Pfaffian version of the classical matrix-tree theorem 
[Kir47j . see also |Cha82] and references therein. Introducing an explicit algorithm, we 
give a constructive proof of our result and a full characterization of half-trees involved. A 
precise statement of our main theorem, as well as consequences for the determinant, are 
given in Section [1.11 of the introduction. An outline of the paper is provided in Section [TT51 
Motivations for proving such a result come from statistical mechanics and are exposed in 
Section 11.31 

1.1 Statement of main result 

Let = V U R, where V = {1, ■ ■ ■ , n}, R = {n + 1, - ■ ■ ,n + r} and n is even. Let 
= {aij)^^ j ev^} be a skew-symmetric whose column sum is zero and let A = {aij)^ij g^} 
be the matrix obtained from A^ by removing the r last lines and columns. The matrix A 
is also skew-symmetric and the Pfaffian of A, denoted Pf(^), is: 

where V is the set of partitions of {1, • • • , n} into n/2 unordered pairs, also known as pair- 
ings, and (Ttt is a permutation such that {(77r(l)(T7r(2), • • • , aTr{n — l)(77r(n)} is a description 
of the partition vr. If n is odd, then by convention Pf(A) = 0. 
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To the matrix one associates the graph = {V^,E^), where = {ij : i,j (z 
V^, Qij 7^ 0}. Every oriented edge (z, j) of is assigned a weight aij, thus defining 
a skew-symmetric weight function on oriented edges. The matrix is the weighted 
adjacency matrix of the graph G^. 

A spanning forest of G^ is an oriented edge configuration of G^, spanning vertices of 
V , such that each connected component is a tree containing exactly one vertex of R. 
This vertex is taken to be the root and edges of the component are oriented towards 
it. Equivalently, a spanning forest of G^ is an oriented edge configuration containing no 
cycles, such that each vertex of V has exactly one outgoing edge of the configuration. A 
leaf of a spanning forest is a vertex with no incoming edge. 

Let G = (y, E) be the graph naturally associated to the matrix A. A perfect matching 
Mq of G is a subset of edges such that each vertex of V is incident to exactly one edge of 
Mq. Note that a perfect matching of G contains exactly \V\/2 edges. In the whole of this 
paper, we suppose that G has at least one perfect matching; if this is not the case, then 
Pf(^) = (see also Section [27T|) . We let M. denote the set of perfect matchings. 

Let F be a spanning forest of G^, then F is said to be compatible with Mq if it consists 
of the \V\/2 edges of Mq and of \V\/2 edges of E^ \ Mq. The oriented edge configuration 
F \ Mq is referred to as a half-spanning forest. In the specific case where R is reduced to 
a point, then i*" is a tree and F \ Mq is referred to as a half-tree. 

Example. Let = {1,2,3,4,5}, V = {1,2,3,4}, R = {5}. Consider the graphs G^ 
and G pictured in Figure [1] below. A choice of perfect matching Mq of G is pictured in 
white, and Fi,F2,F^ are examples of spanning trees of G^ compatible with Mq. 




F, F, 
Figure 1: Spanning trees compatible with Mq. 

Here is the statement of our main theorem. 

Theorem 1 (Principal minors Pfaffian half-tree theorem). Let A^ be a skew- symmetric 
matrix of size (n + r) x (n + r), whose column sum is zero, such that n is even; and let A 
be the matrix obtained from A^ by removing the r last lines and columns. Let G^ and G 
be the graphs naturally constructed from the matrices A^ and A, respectively. 

For every perfect matching Mq of G, the Pfaffian of A is equal to: 

Pf(^) = ^ sgn{aMo{F\Mo)) n 

where Og is the coefficient of the matrix A^ corresponding to the oriented edge e; sgn{an.fg(^p\n.fg^) 
is the signature of the permutation cr]\,f^^^p\]\,f^^-j of Definition [7771 below; F{Mq) is the set 
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of spanning forests of compatible with Mq, satisfying Condition (C) of Definition M.^ 
below. 

Definition 1.1. Let F be a spanning forest of compatible with Mq. The orien- 
tation of F induces an orientation of edges of the perfect matching Mq, and we let 
(^1)^2), ■ ■ ■ , (in-ijin) be a description of the oriented matching. Then, o"jvfQ(_F\Afg) is the 
permutation: 

_ f 1 2 ••• n 

Note that the interchange of two pairs does not change the signature of the permutation. 

Here is the algorithm used to characterize half-spanning forests which contribute to Pf(^). 
Trimming algorithm 

Input: a spanning forest F of G^ compatible with Mq. 
Initialization : Fi = F. 
Step i, i > 1 

Let £\ be the largest leaf of Fi and consider the connected component containing i\. Start 
from £\ along the unique path joining i\ to the root of this component, until the first time 
one of the following vertices is reached: 

- the root vertex, 

- a fork, that is a vertex with more than one incoming edge, 

- a vertex which is smaller that the leaf i\. 

This yields a loopless path A^i starting from £\, of length > 1. Let Fj+i = Fi\ A^i . If Fj+i 
is empty, stop; else go to Step i + 1. 

End : since edges are removed at every step, and since F contains a finite number of edges, 
the trimming algorithm ends in finite time A^. 

Definition 1.2. A spanning forest F compatible with Mq is said to satisfy Condition (C) 
if each of the paths A^i , • • • , A^jv obtained from the trimming algorithm, starts from an 
edge of Mq and has even length. We let J^{Mq) be the set of spanning forests compatible 
with Mq, satisfying Condition (C). 

Example. Applying the trimming algorithm to each of the spanning forests i^i, F2, F3 of 
Figured] yields: 



Fi 
F2 
F3 



Step 1 
Step 1 
Step 1 



?i = 2,A2 = 2,3, L Step 2 
?i = 2, A2 = 2,3,5. Step 2 
f\ = 4, A4 = 4, 1. Step 2 



ff2 



f = l,Ai = 1,4,5. 
il = l,Ai = 1,4,5. 
tj = l,Ai = 1,2,3,5. 



The spanning trees -Fi and F2 satisfy Condition (C) but not F-^ 



3- 
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Remark 2. 



• It is interesting to note that taking different perfect matchings Mq yields different 
families of half-spanning forests. It is not clear a priori, without using the Pfaffian 
half-tree theorem, that these families should have the same total weight. 

• Suppose that we change the labeling of the vertices. Let be the corresponding 
skew-symmetric adjacency matrix, and A be the matrix obtained by removing the r 
last lines and columns. As long as the re-labeling does not affect vertices of R, the 
matrix A is obtained from the matrix A by exchanging lines and columns, so that 
Pf(A) = ±Pf(^). Applying the Pfaffian half-tree theorem to the matrices A and A 
nevertheless yields a different set of half-spanning forests and again, it is not clear 
a priori that they should have the same total weight in absolute value. Note that 
taking other principal minors amounts to changing the labeling of the vertices. 

• In the paper [MV02| . Masbaum and Vaintrob assign to a 3-uniform hypergraph 
a specific skew-symmetric matrix whose column sum is zero, and prove that the 
Pfaffian of any principal minor of this matrix enumerates signed spanning trees of 
the 3-uniform graph. A combinatorial proof of their result is then given by Hirschman 
and Reiner |HR04] , and yet another proof using Grassman variables is provided by 
Abdesselam |Abd04) . The matrix considered by Masbaum and Vaintrob satisfies the 
assumptions of Theorem [H and it would interesting to investigate whether there is 
an explicit relation between half-spanning forests of Theorem [1] and spanning trees 
of the 3-uniform graph. 

Using the fact that the determinant of a skew-symmetric matrix is the square of the 
Pfaffian, we obtain the following corollary. 

Corollary 3. Let A^ be a skew-symmetric matrix of size (n + r) x {n-\-r), whose column 
sum is zero, such that n is even; and let A he the matrix obtained from A^ by removing 
the r last lines and columns. Let and G be the graphs naturally constructed from the 
matrices A^ and A respectively. 

The determinant of the matrix A is equal to: 

det(A)= Yi E n«-' 

Mo ex FeT{Mo) eeF 

where Oe is the coefficient of the matrix A^ corresponding to the oriented edge e, and 
F{Mq) is the set of spanning forests compatible with Mq, satisfying Condition (C). 

Remark 4. 

• The fact that principal minors of a skew-symmetric matrix whose column sum is 
zero, count spanning forests is also a consequence of the more general all-minors 
matrix-tree theorem (which holds for any matrix whose column sum is zero). A 
combinatorial way of proving this result is to use the explicit expansion of configura- 
tions due to Chaiken |Cha82] . This method is not satisfactory in our context, since 
it does not shed a light on how spanning forests are obtained from double perfect 
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matchings, which is what we aim for, see Section [1.3[ Indeed, the idea of Chaiken's 
proof is to expand terms on the diagonal of the matrix and show that only spanning 
forests remain. In the case of skew-symmetric matrices, since diagonal terms are 
this amounts to 'artificially' creating configurations which do not exist. As a result 
of our proof, we explicitly construct spanning forests from double perfect matchings, 
and identify a specific family of spanning forests counted by principal minors. In 
particular, this implies that in the case of skew-symmetric matrices, specific cancel- 
lations occur within spanning forests of the general matrix-tree theorem, a fact hard 
to establish without using Corollary [3l 

• An intrinsic definition of UA/ggyn J-'(Mo), not using reference perfect matchings, is 
given in Remark [16] of Section 13. 2i 

• A line bundle version of this result, in the spirit of [For93] and [Kenllj . is proved in 
Section [331 see Corollarv [T9l 

1.2 Outline of the paper 

• Section [21 In Section 12.11 we state the interpretation of the Pfaffian as counting 
signed perfect matchings of the graph G. Fixing a reference perfect matching Mq, 
we then introduce an explicit algorithm, which constructs from the superimposition 
of Mq and a generic perfect matching M counted by the Pfaffian, a family of half 
RC-spanning forests whose connected components are trees rooted on vertices of R, 
or on cycles of even length > 4; and whose total weight is equal to the contribution 
of M to the Pfaffian. The main tool of the algorithm is the 'opening' of doubled 
edges procedure, described in Section 12.31 Step 1 of the algorithm is exposed in 
Section [2.41 and the full algorithm is the subject of Section [2.5i A characterization 
of configurations obtained is given in Section 12. 6[ 

• Section 3 . Section 13.11 consists in the proof of Theorem [H The idea is to show that 
the contribution of half RC-spanning forests constructed above, having connected 
components rooted on cycles of length > 4 cancel, and that only the contribution of 
spanning forests (rooted on vertices of R) remains. The characterization of configu- 
rations obtained from the algorithm is also simplified in the case of spanning forests, 
yielding the trimming algorithm of Section 11.11 of the introduction. The proof of 
Corollary [3l is the subject of Section [3^21 Finally, in Section [3^ Corollarv [T9l proves 
a line bundle version of the matrix-tree theorem for skew-symmetric matrices of 
Corollary [3l 

1.3 A question from statistical mechanics 

The Pfaffian half-tree theorem answers a question raised when working on the paper |dT] , 
in which we prove an explicit relation, on the level of configurations, between two mod- 
els of statistical mechanics: the dimer model on the Fisher graph corresponding to the 
low temperature expansion of the critical Ising model (through Fisher's correspondence 
|Fis66] ). and spanning forests. The question raised does not rely on the Fisher graph and 
can be rephrased in the following, more general framework. 
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In the setting of statistical mechanics, a perfect matching of a graph is known as a dimer 
configuration. Assigning non-negative weights to edges of the graph naturally defines a 
weight for each dimer configuration (by taking the product of the edge-weights present 
in the configuration) and a probability measure on all dimer configurations of G, thus 
yielding a statistical mechanics model. The dimer model on planar graphs has been the 
subject of extensive studies in the last 50 years, and of huge progresses in the last 15 years, 
see |Ken09] for an overview. 

A double dimer configuration is the superimposition of two dimer configurations. Our 
goal is to explicitly construct spanning forests from double dimer configurations when the 
model is critical, and to do so on the same graph, thus relating two models of statistical 
mechanics on the level of configurations. 

When the graph is planar, dimer configurations are counted by the Pfaffian of the Kaste- 
leyn matrix |Kas67t ITFGlj . which is a skew-symmetric matrix with a specific choice of 
orientation of the edges compensating signs of signatures. The square of the Pfaffian, 
which is the determinant, counts double dimer configurations. In the case of the dimer 
model corresponding to the critical Ising model the column sum is zero (when the graph 
is embedded on the torus), a fact related to the model being critical. This motivates our 
choice of taking the column sum equal to zero. In general, the Pfaffian of a skew-symmetric 
matrix counts signed perfect matchings, with the sign coming from the signature of a per- 
mutation naturally assigned to the matching, see Section 12.11 Theorem [1] proves that 
principal minors of the Pfaffian of a general skew-symmetric whose column sum is zero 
counts a specific family of half-spanning forests, and half-spanning forests are explicitly 
constructed from perfect matchings. Corollary [3] proves that principal minors of the deter- 
minant of such a matrix counts a family of spanning forests, and the latter are explicitly 
constructed from double perfect matchings. Specifying this result to the case of planar 
graphs or graphs embedded on the torus, and Kasteleyn matrices, answers our initial 
question. 

It is known from the work of Temperley [Tem72j . Kenyon, Propp and Wilson [KPWOO] 
that spanning trees are in bijection with dimer configurations of a related bipartite graph. 
Conversely, Kenyon and Sheffield |KS04j relate dimer configurations of a bipartite graph 
to spanning trees of a process on an associated graph. In this paper perfect matchings and 
trees live on the same graph, the graph must not be bipartite nor even planar, the weight 
function on edges of the graph must not be positive, but the column sum must be zero. 
The main drawback of our result in the context of statistical mechanics is that, even when 
the matrix is Kasteleyn and perfect matchings all have positive weights, corresponding 
spanning forests might have negative weights. 

2 Prom matchings to half i?C-rooted spanning forests 

Let us recall the setting: is a skew-symmetric matrix of size {n + r) x (n + r), whose 
column sum is zero, such that n is even; and A is the matrix obtained from A^ by removing 
the r last lines and columns; and G are the graphs naturally constructed from the 
matrices A^ and A in Section 11.11 of the introduction. 
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Definition 2.1. An RC -rooted spanning forest, referred to as an RCRSF is an oriented 
edge configuration of spanning vertices of G, such that each connected component is, 
either a tree rooted on a vertex of R, or a tree rooted on a cycle of G, which we refer to 
as a unicycle. Edges of each of the components are oriented towards its root, and edges 
of the cycles are oriented in one of the two possible directions. 

Definition 2.2. Let Mq be a reference perfect matching of G. An RCRSF F is said to be 
compatible with Mq, if it consists of the \ V\/2 edges of Mq, and of \V\/2 edges of \ Mq. 
Moreover cycles of uni-cycles have even length > 4, and alternate between edges of Mq 
and F \ Mq. The oriented edge configuration F \ Mq is referred to as a half-RCRSF. 

In Section 12.11 we give the graphical interpretation of the Pfafhan of the matrix A as 
counting signed perfect matchings of G. Let M be a generic perfect matching counted by 
the Pfaffian and Mq be a fixed reference perfect matching of G. In Sections 12.41 and \2.5\ 
we introduce an explicit algorithm which constructs, from the superimposition of Mq and 
M, a family of half i?C-rooted spanning forests compatible with Mq, whose total weight 
is equal to the contribution of M to the Pfafhan. In Section 12. 6| we characterize RG- 
spanning forests obtained. Notations used are given in Section 12.21 The main graphical 
idea of the algorithm is the subject of Section [231 

2.1 Graphical interpretation of the PfafRan 

Recall that V denotes the set of pairings of {1, • • • ,n}, and let M be the set of perfect 
matchings of G. Observing that every perfect matching of G corresponds to a pairing of 
V, and that pairings of V which do not correspond to perfect matchings of G contribute 
to the Pfaffian, we can rewrite Pf(^) as: 

Pf(^) = Sgn(fTM)a<7M{l)fTM{2) • • • a<7M{n-l)'TM{n)' 

MeM 

where au is a permutation such that {(Tm (1)cm(2), • • • , (yM{n — l)c'"Af (f^)} is a description 
of the perfect matching M. 

Choosing the permutation gm amounts to choosing an order for the n/2 pairs and an 
order for the two elements of each of the pairs, meaning that there are (^)!22 choices. 
Exchanging two pairs does not change the signature nor the corresponding coefficients 
of the matrix, whereas changing two elements of a pair changes the signature of the 
permutation and the sign of the corresponding element of the matrix. As a consequence, 
the global sign is unchanged, and fixing the signature amounts to choosing an orientation 
of edges of the perfect matching. 

We now specify the choice of signature of the permutation cjm by choosing an orientation 
of edges of M. Let Mq be a fixed reference matching of G. The superimposition of Mq 
and M, denoted by Mq U M, consists of disjoint doubled edges (covered by both Mq 
and M) and alternating cycles of even length > 4, covering all vertices of G. Let us, 
for the moment, consider doubled edges as cycles of length 2. The orientation of the 
superimposition Mq U M is fixed by the following rule: for each cycle, the orientation is 
compatible with that of the edge where is the smallest vertex of the cycle, and 
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i'l is its partner in Mq. This yields an orientation of edges of M, and thus a choice of cjm- 
This procedure also gives an orientation of edges of Mq, and thus a choice of (Tmo- Let us 
also denote by Mq U M the oriented superimposition. 

Example (Figure [2]). In white is a choice of reference perfect matching Mq and in black 
are the three possible matchings Mi, M2, M3 of the graph G. Edges of the respective 
superimpositions are oriented according to the rule described above. 

M :q; 

'■■•■5 '■•■■5 '■•■■5 

Ml Mj M3 

Figure 2: Oriented superimposition. 



Let (TMoUA/ be the permutation whose cyclic decomposition corresponds to cycles of the 
superimposition Mq U M. Then, 

Sgn(aMoUM) = (-l)l^(^oUM)|(_i)|C{MoUM)|^ 

where V{Mo U M) is the set of doubled edges of Mq U M and C{Mq U M) is the set of 
alternating cycles of length > 4 of Mq U M. Note that this signature does not depend on 
the orientation of the cycles. Following Kasteleyn |Kas67| . signatures of the permutations 
(TMo and (Tm are related as follows: 

sgn((TM) = sgn(fJMo)sgn(fJMoUM) 

= sgn(aMo) (-l)l^(A^oUM)|(_i)|C(MoUM)|_ 

Writing aMo as to remember that our choice of orientation of edges of Mq depends 

on M, the Pfafhan of A can be expressed as: 

Pf(A)= ^ wmAM), (1) 
MeM 

where 

r.Mo(M) =sgn(aMo(M))(-l)l^(^'°^''')l(-l)l"(''°^''')l H 

egM 

and Oe is the coefficient of the matrix A corresponding to the oriented edge e of M. 



2.2 Notations 

Let Mq be a fixed reference perfect matching of the graph G, M be a generic perfect 
matching, and Mq U M be the oriented superimposition of Mq and M constructed in 
Section 12.11 In order to shorten notations, we write T> instead of 'D{Mq U M) for the set 
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of doubled edges of the superimposition, C instead of C{Mq U M) for the set of cycles of 
length > 4, and w{M) instead of wmo{^)- 

We now introduce definitions and notations used in the algorithm of Sections 12.41 and 12.51 
Let denote the set of vertices of V which belong to a cycle of C. For every subset T>' of 
doubled edges of V, let denote the set of vertices of V which belong to doubled edges 
of v. 

Every vertex i G belongs to a doubled edge covering vertices i and i' oiV. We denote 
by ej (or this doubled edge and define to be the set of vertices in the full graph G^, 
adjacent to i' other than i. 

For every subset V ofV, denote by the set of vertices of Vi which belong to doubled 
edges of D', and by (V^ Y those which don't. Then Vi can be partitioned as: Vi = 

vr u {vr'r. 

2.3 Idea of the algorithm 

The idea of the algorithm is to use the reference configuration Mq as a skeleton for opening 
up doubled edges of the superimposition Mq U M. Indeed, because of the condition 
X]j=i — 0' configurations of Figure [3] have opposite weights. 




Figure 3: 'Opening' of doubled edges procedure: ai' i = — ^ Oj' j. 

There are two main difficulties in realizing this procedure: the first is that there is, a 
priori, no natural way of deciding whether to 'open' up of the doubled edge at the vertex i 
or at the vertex i' . The second is that we want to keep track of configurations constructed, 
show that we obtain i?C-rooted spanning forests, characterize them and prove that only 
spanning forests remain. It turns out that the 'opening' procedure depends strongly on 
the labeling of vertices. 

2.4 Algorithm: Step 1 

Recall that the goal of the algorithm is to construct, from the superimposition Mq U M 
of a reference perfect matching Mq and a generic perfect matching M, a family of half- 
i?C-rooted spanning forests of compatible with Mq, whose total weight is equal to 
the contribution of M to the Pfaffian. In this section, we introduce the first step of the 
algorithm, setting rules for the opening up of doubled edges of MqUM. The full algorithm, 
which in essence consists of iterations of Step 1, is the subject of Section [2.51 

Input: oriented superimposition Mq U M. 
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Initialization : if the superimposition Mq U M consists of cycles only, that is if the set T> is 
empty, let Oq = {M} and stop. Else, let Oq = {0} and go to the first iteration. 

Example (Figure d]). Consider Figure [3] as input of the algorithm. The algorithm will 
be explicitly performed on this example, throughout Sections 12.41 and | 2.5[ 




'■■■■5 



Figure 4: Input Mq U M of Step 1. 

Since Mq U M contains doubled edges, the output Oq is {0}. 
Iteration 1 

• Define li = min{z E V : i belongs to a doubled edge of P}. Then ii is the partner 
of a vertex i'^ in Mq and M. By our choice of orientation for Mq and M, the edge 
iii'i is oriented from ii to i'l in Mq and from i[ to ii in M. For every £2 E V^^, 
define: 

Me„e, = {M\{i[,ei)}u{{e[,£2)} 
t/;(M,,,,,) = sgn(c7Mo(Af))(-l)"''"'(-l)"'' J] 

Example (Figure (5]): ^1 = 1, Ve^ = Vi = {2,3,5}, yielding configurations Mi_2, 

Mi,3, Mi,5. 




Figure 5: From left to right: black edges are the oriented edge configurations Mi^2, -^1,3) 

Mi,5. 



• Let be the set of doubled edges V \ {e^^}. Then, the set V^^ can be partitioned 
as the set of vertices of V^j which belong to a doubled edge of Vi^ and the set 
of those which don't. Using notations of Section 12.21 this can be rewritten as: 

y,, = y,^u(y,^r. 
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The output of Iteration 1 is the set of configurations M^^^^i^ such that £2 does not 
belong to a doubled edge of T>£-^ : 

w{0{) = Yl 

where by convention, if (V^^ = 0, then Oi = and w{Oi) = 0. 

• The algorithm continues with configurations Mi^^i^ where £2 belongs to a doubled 
edge of P^j. Formally we have: if V^^ ^ = 0, then stop; else, go to Iteration 2. 

Example: the set = Vi consists of the doubled edge 23. As a consequence, 
the set Vi-^ = Vi = {2,3,5} is partitioned as Vi = {2,3} U {5}, and the output of 
iteration 1 is Oi = {Afi^s}. The algorithm continues with Mi^2 and Mi^^. 

Iteration k, {k > 2) 

For every £2 G V^'^ , • • • ,£k e V^^^^'" ''^'"^ , do the following. 

• The vertex is the partner of a vertex in Mq and M (since P^^^... is a 
subset of V). If £k < £'/., then by our choice of orientation, the edge £k£'i^ is oriented 
from £k to in Mq and from to £k in M^^^...^^^,. If £fc > then we change 
the orientation of this edge in Mq and in M^^ Let us also denote by Mi-^ ...^£^ 
this new configuration. This change of orientation has the effect of changing the 
permutation assigned to Mq, and we denote by ^ j this new permutation. 
It also negates the contribution of M^j^...^^^, so that the global contribution remains 
unchanged. For every ^fc+i G V^^ , define: 

Me,,..,i,^, = (M,,,..,4 \ (4,4)) U (4, 4+1) 
«;(M,,,..,,,^J = sgn(aMo(Af,,,..,j)(-l)l''l-'(-l)l''' H (2) 

Example (Figure [6]). Recall that V^'^ = V^'^^^ = {2,3}, so that 4 G {2,3}. If 
4 = 2, then = 1^2 = {1,4,5}, yielding configurations Mi^2,i, ^1,2,4, ^1,2,5- If 
4 = 3, then = V3 = {1,4}, yielding configurations Mi^s^i, Mi^3^4. 

• Let ... be the set of doubled edges \{e^^}. Then, the set Vg^ can be 

partitioned as: Yi^ = V^^ ^ U(^;. ^' '°)'^' and the output of Iteration k is the 
set of configurations M^^^... ^^^.^^ such that ^fc+i does not belong to a doubled edge of 
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'■■'5 '■■■'5 



Figure 6: First line, from left to right, black edges consists of the configurations Mi^2,i, 
^1,2,4, ^1,2,5- Second line, from left to right, black edges consists of the configurations 



Ok= [J ■■■ U U ^^i-",4+i, 

• If Vg^^^'"''"'' = 0, then stop. Else, go to Step k + 1. 

Example: when £2 = 2, the set T>i^2 is empty so that V2 is partitioned as {0} U 
{1, 4, 5} and the contribution to the output O2 of iteration 2 is Mi^2,i, -^i,2,4j ^^1,2.5- 
When £2 = 3, the set Pi^s is also empty, implying that V3 is partitioned as {0}U{1, 4} 
and the contribution to the output O2 of iteration 2 is Mi^s,i, -^1,3,4- After Iteration 
2, for every £2 € V^^^S the set V^^^^^^ empty, so that the algorithm stops. 

End 

Step 1 of the algorithm stops at time m for the first time, if it hasn't stopped at time 
m - 1, and if for every £2 € vl^'\- ■ ■ .£„, € y/^'i'- 1 ; v?'^'"'''"' = 0. This implies in 

particular that (V;,^'' ' ' = V^^. Since the number of doubled edges decreases by 1 
every time an iteration of the algorithm occurs, and since the number of doubled edges in 
V is finite, we are sure that Step 1 of the algorithm stops in finite time. 

2.4.1 Output of Step 1, geometric properties of configurations 

The output of Step 1 of the algorithm is the set of configurations Si = UfcLo^fc- The 
weight of this set is defined to be w{Si) = Ylk=o''^i^k)- 
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If the initial superimposition Mq U M consists of cycles only, i.e. if the set D is empty, 
then m = and Si = {M}. In all other cases, the set Si can be rewritten in a more 
compact way as: 

where: 

r^i = {7^1 : lii is a path of length 2k for some A; e {1, • • • , m} : 7^, = h^l'i, ■ ■ ■ ,4,4)4+1, 
li = min{i : i belongs to a doubled edge of D}, 



Vj G {2, • • • , k], £j G V^_\' ' ' and is the partner of ij in Mq and M, 



Let 7£j = - ,^fc,4,4+i be a generic path of F^^ for some k G {!,••• ,m}, and 

let -F^^^ denote the superimposition Mq U M^^^ . The configuration F^^^ and the path 7^1 
satisfy the following properties. 

• The oriented edge configuration F^^^ : 

(I) has one outgoing edge at every vertex of V, and contains the path 7^^. 
(II) has k doubled edges less than Mq U M. 

• The oriented path 7^^ : 

(III) has even length 2k, is alternating (meaning that edges alternate between Mq and 
M-y^^). It starts from the vertex £1 followed by an edge of Mq. 

(IV) The vertex ii is the smallest vertex belonging to a doubled edge of D. The 2k 

first vertices of 7^^ are all distinct and the last vertex 4+i belongs to (V^^ ^' ' *)^. 
Observing that: 

K'""''''r = Vi, n (i?u y'^ u {£i,e'i, ■ ■ ■ ,4,41), 

we deduce that one of the following holds. 

(IV) (1) If 4+1 G R, then 7^^ is a loopless oriented path from £1 to one of the root vertices 
of R, and £1 is a leaf of F^^^ . Since R = {n + 1, ■ ■ ■ ,n + r}, £1 is smaller than all 
vertices of 7^^ . 

(IV) (2) If £k^i G V^, then 7^^ is a loopless oriented path ending at a vertex of one of the 
cycles of C that is, the connected component containing £1 is a unicycle with a unique 
branch. The vertex £1 is a leaf of F^^^ and is smaller than the 2k first vertices of the 
path, but cannot be a priori compared to vertices of the cycle of C By construction of 
the orientation of Mq U M, see Section [2.11 the orientation of the cycle is compatible 
with that of the edge {11,12), where ii is the smallest vertex of the cycle and ^2 is its 
partner in Mq. 
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(IV)(3) If ik+i € {ii,£'i, ■ ■ ■ then 7^^ contains a loop of length > 3. If i^+i = U for 

some i € {!,••• , A;}, then the loop has even length and is alternating and the part 
of to the loop also has even length, is alternating and starts with an edge of Mq. 
Moreover, the orientation of the loop is compatible with the orientation of the edge 
the vertex ^1 is smaller than all vertices of the path to the cycle and 
smaller than all vertices of the cycle. 

Note that if Ik+i 7^ ^1, then £1 is a leaf and the connected component containing (.1 
is a unicycle with a unique branch. Else, if H-k+i = ^1, the connected component is a 
cycle. 

If ik+i = (■[ for some i G {!,••• , A:}, then the loop has odd length with two edges of 
M incident to the vertex I'-. Observing that the loop in both directions is obtained 
from Step 1 of the algorithm, and using the fact that the matrix A is skew-symmetric, 
we deduce that the contributions of these configurations cancel and we remove them 
from the output of Step 1. Thus we only consider configurations such that ^a:+i = f-i 
for some i € {1, • • • , /c}. 

Example. The output Si of Step 1 is: Si = {Mi,5, Afi,2,i, Mi,2,4, Afi,2,5, Mi,3,i, Mi,3,4}. 
Configurations Mi^s, Mi^2,5 are in Case (IV)(1), configurations Mi^2,i, Afi,3,i are in Case 
(IV) (3) with even cycles created, and configurations Mi^2,4; -^1,3,4 are in Case (IV) (3) with 
odd cycles created. Contributions of Mi ^2,4 and Mi ^3^4 cancel so that they are removed 
from the output. As a consequence the final output of Step 1 is, see also Figure [T) 

Si = {Mi,5,Mi,2,l,Mi,2,5,Mi,3,i}, 




Figure 7: Output of Step 1 of the algorithm. 



2.4.2 Weight of configurations 

As a consequence of the next two lemmas, we obtain that Step 1 of the algorithm is weight 
preserving i.e. w{Si) = w{M), see Corollary [71 

Lemma 5. We have: 

. w{M) = HM,,^i2)- 
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• If inn > 2, then for every k € {2, • • • , m} and every £2 G ; ■ ' ' ; ^fc € V^^ ^ : 

^M,,,..,,,) = ^ u;(M,,,..,4^J. 

Proof. Suppose that m >2, the proof in the other case being similar. For every ik+i G Vi^ , 
M,,,..,4^, = {M,,,..,4 \ (4,4)} U {(4, 4+1)}, thus: 

ae = — Be- 

Ofi f, 

By assumption, coefficients of each hne of the matrix sum to 0. Returning to the 
definition of Vf., this imphes that > api = —ap' e,. Thus, 

E n ae=- n (3) 

Combining Equation ^ with the definition of the weight of configurations given in Equa- 
tion yields: 

u;(M,,,..,,,^J = sgn(aM„(M,,,...,j)(-l)l'^l(-l)l^l-^ E 11 «^ 

= sgn(aM„(M,,,...,j)(-l)l'^l(-l)l^|-('=-^) n «^ 

□ 

Lemma 6. Suppose m > 2. Then for every k € {2, • • • , m}, 

m 



i=k 



Proof. In order to simplify notations let us write, only in this proof, instead of 
Vf ^' ' Lemma [6] is proved by backward induction on k. 

Suppose k = m. By definition of the last step of the algorithm, Vj^ = 0, so that {Vj^Y — 
Vl^ and: 

w{Om)= J2 E E 

1 1 t™ — 1 



Y ■■■ E (by Lemma ED, 



^1 *^m— 1 
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thus proving the case k = ni. Suppose that the statement is true for some k € {3, • • • , m}. 
By Iteration A: — 1 of Step 1 of the algorithm, we know that: 



Combining this with the induction hypothesis yields: 

m m 
i=k—l i=k 

= E ••• E ( E + E 

= E E (by LemmaO, 

proving the statement for k — 1 and ending the proof of Lemma El □ 
Corollary 7. 

w{Si) = w{M). 

Proof. Suppose m > 1. Then, 

m 
k=2 

= w{Oi)+ Y w{Me,/,), (byLemmaED 

= Yl wiMi,^e^)+ Y ^(^^a), (by definition of Oi) 

= Y ^(^4/2) 

= w{M), (by Lemma [5]). 
When m = 0, Si = {M}, and the conclusion is immediate. □ 

2.5 Full algorithm 

Let Mq be a reference perfect matching of the graph G and let M be a generic one. Recall 
that C denotes the set of cycles of length > 4 of the superimposition Mq U M, and T> 
denotes its set of doubled edges. 
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In Section 12.41 we established Step 1 of the algorithm, starting from the superimposition 
Mq U M, yielding a set of oriented edge configurations Si through the opening of doubled 
edges procedure, whose total weight is equal to the contribution of M to the Pfaffian. 
In this section, we introduce the full algorithm, which in essence consists of iterations of 
Step 1 performed until no doubled edges of V remain. 

Let us directly handle the following trivial case. If Mq U M consists of cycles only, that is 
if the set V is the empty, then the opening of edges procedure does not start, and recall 
that the output of Step 1 is 5i = {M}. The same holds for the full algorithm and its 
output is r = {M}. 

2.5.1 Step 1 of the full algorithm 

Assume that the initial superimposition contains at least one doubled edge, i.e. T> ^ 0. 

Notations are complicated by the fact that the algorithm depends on the labeling of the 
vertices and that iterations of Step 1 depend on previous steps. We thus need many indices 
to keep track of everything rigorously, but one should keep in mind that, in essence, we 
are iterating Step 1. Let us add sub/superscripts to Step 1 of Section [231 That is l\ 
becomes £}, Iteration k becomes k\ and Step 1 ends at time m\. The set of configurations 
obtained from Step 1 is, S\ = \J ^ ^ifi ' weight is w{Si) = ^ ^ w(My^^ ). 

For every 7^1 € L^i , let P^^^ be the set of doubled edges of the superimposition Mq U M^^-^ . 
If 1=0, then stop; else go to Step 2. 

Output of Step 1 of the full algorithm. It is the subset 7i of 5i, consisting of configura- 
tions -A^7^i where 7^1 € F^i, and P^^j is empty. Formally, 

ri= U M^.r 

7fi6(r^i)<= 

where F^i is the set of paths 7^1 of F^i such that P^^^ is non-empty. If for all 7^1 € F^i, 
the set ^ is non empty, then 71 = 0. 

Example. Recall that the output of Step 1 is 5i = {Mi^s, Mi^2,i; ^1,2,5, ^1,3,1}- The 
set of doubled edges of the superimposition of Mq and Mi^2,i; ^1,2,5, -^1,3,1 is empty, so 
that the output of Step 1 of the full algorithm is, 71 = {Mi_2,i, Mi^2,5, -^i,3,i}i and the 
algorithm continues with the configuration Mi_5. 

2.5.2 Step j of the full algorithm, j > 2 

For every 7^1 G r«i,--- ,7„j-i G F.j-i (7^1 , • • • ,^.3-2), perform Step 1 of the algorithm 

with the initial superimposition Mq U . That is, define i\ = min{z E V : 

1 *i 

i belongs to a doubled edge of D^^^^... }, and iterate until the algorithm ends at some 
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time rrij. Everything works out in the same way because P^^^^... ^-y is a subset of T>. 

1 ^1 

The output is the set of oriented edge configurations: 
where 

^4(7^1, • • • >7£i-i) = ^lei ■ ^ei is ^ P^*^ ^i'-^i ' ' ' ' '^i'^i'^i+i' ^^me fe^ G {1, • • • ,mj}, 

such that = min{z EV:i belongs to a doubled edge of Vj i,...,^ j.^} 

1 ^1 

and ij is the partner of in Mq and M, 

V i i \ 

7,l,-,7„.,-l/i,-.4,. 



^1+1 ^ 



For every G r^j(7^i, • • • do the following: if V^^^^... is empty, stop; else go 

to step J + 1. 

Output of Step j of the; I'ull algorithm. 

Let 7} (7^1, ■ ■ ■ )7£j-i) be the subset of 5,(7^1, • • • , T^j-Oj consisting of configurations My^^^...^^^. 
such that ■ is empty. Formally, 

1 ^1 

"T^- (74 , • • • , 7^.^-1 ) = ^ U ,7,.^ , 

where F^j (7^1 , • • • , 7^3-1 ) is the subset of paths 7^j of F^j (7^1 , • • • , 7£i-i ) such that "D-y^j 
is non-empty. If X>-y^j^...^^^^ is non empty, then 7j(7^i, • • • )7£j-i) = 0- 

Then, the output Tj of Step j of the full algorithm is: 

7^= U U U ^7,J,-,7,,- 

For convenience, we shall also use the notation Fj for the paths (7^1 , ■ ■ ■ , 7^^ ) involved in 
Tj, i.e: 

rj= U ^74,-,74- 

(7<.i,->7^j)erj ^ 



18 



The weight of Tj is the sum of the weights of the configurations it contains. 

Example. In Step 2, we perform Step 1 of the algorithm starting from the initial su- 
perimposition Mo U Mi^5. The latter contains one doubled edge 23, thus the vertex 
i\ is the smallest of 2 and 3, that is 2. The output 52 consists of the configurations 
-^i,5;2,i, -M'i,5;2,4, Mi^5;2,5, depicted in Figure [8] below. 




Figure 8: Output of Step 2 of the algorithm. 

The superimposition of Mq and the above three configurations contains no doubled edges. 
As a consequence, the full algorithm stops and the output T2 of Step 2 is: 

T2 = {Mi,5;2,l,Mi,5;2,4,Mi,5;2,5}. 

2.5.3 End and output of the full algorithm 

The algorithm stops at Step T for the first time, if it hasn't stopped at time T — 1, and 
if for every 7^1 G F^i,--- ,7^t-i G F^t-i (7^1 , • • • ,1^-2) , 7^t G F^t(7£1,--- the 

11 ^1 11 1 111 ^1 

superimposition MoUM-y i,...,^ j. contains no doubled edge. This implies in particular that 

{TfT (7^1 , • • • , lf-i)Y = TfT (7^1 , • • • , 7^^-! ) . Since the number of doubled edges decreases 
at every step and since T) is finite, we are sure that this happens in finite time. 

The output T of the full algorithm is : 

T 

J = l 

Example. The output of the fuh algorithm is: 

r =riUr2 = {Mi,2,1,Mi,2,5,Mi,3,i,Mi,5;2,1,Mi,5;2,4,Mi,5;2,5}, 

summarized in Figure O below. 

The weight of T is the sum of the weights of the configurations it contains. If the initial 
superimposition Mq U M consists of cycles only, i.e. the set T> is empty, then T = {M}, 
and 

w{T) = w{M) = sgn(aMo{M))(-l)'^' H 

eGM 
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M, 




3 4 




2 1 



3 4 



2 1 



3 4 




M, 



M, 



M, 



M, 



Figure 9: Output of the full algorithm. 



In all other cases: 

T 

i=i (74,-,7,j)6rj 
1 '^i 

Since for every j G {!,••• j^}; and for every (7^1,- '' jT^^) ^ ^j, the superimposition 

Mq U My J. contains no doubled edge of V, we have: 

1 *i 

"'(^7,l,-,7,.) = Sgn(fTMo(M^^^,..,^^.))(-l)''^' n 

1 «i 

By iterating the argument of Section 12.4.21 we obtain the following: 
Corollary 8. The full algorithm is weight preserving, that is: 

w{T) = w{M). 

2.6 Geometric characterization of configurations 

Consider the superimposition Mq U M, recall that C denotes the set of cycles of length 
> 4 of Mq U M, and that V denotes its set of doubled edges. Consider the full algorithm 
with initial superimposition Mq U M in the case where Mq U M contains doubled edges, 
that is when D ^ 0. Let j £ {1, • • • (7^1 , • • • ,7^j) G Tj, and M^ i,- -,7 j G T be a 

generic output; and denote by ^ . the superimposition Mq U My ^ . In order to 

simplify notations, we introduce: 

Vi e {!,••• F,:=F ... 

One should keep in mind that the index j refers to the last step of the algorithm, and that 
indices i G {1, • • • ,j — 1} refer to intermediate steps. As a consequence of the algorithm, 
see (HD, the configuration Fj and the paths 7^1, •• • ,7^j satisfy the following properties. 

• The oriented edge configuration Fj: 

(I) has one outgoing edge at every vertex of V. It consists of the paths 7*1 , • • • , "Jp and 

1 t-i 

of the cycles C of the initial superimposition Mq U M; 

(II) has no doubled edge of V since the full algorithm precisely stops when this is the 
case. 
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For every i £ {!,■ ■ ■ the path satisfies the following. 

(III) It has even length 2ki for some ki E {1, • " " j^t-j} and is alternating. It starts from 
the vertex £\ , followed by an edge of Mq . 

(IV) The vertex £\ is the smallest vertex belonging to a doubled edge of (understood 
as Mq U M when i = 1). The 2ki first vertices are all distinct and the last vertex 
'^ki+i belongs to: 

As a consequence, one of the following 5 cases hold. 

• If e -R U U {4, (fi)', • • • , 4i > (^L)'U ' 7^1 consists of a new connected 
component, and we recover the three cases obtained after Step 1 of the algorithm, 
replacing 7^^ by , see Section [2.4.1l For convenience of the reader, we repeat these 
cases here. 

(IV) (1) If 4._|_i S then 7^2 is a loopless oriented path from £\ to one of the root 
vertices of R. Since R = {n + 1, • • • , n + r}, l\ is smaller than all vertices of 
7^i . The vertex i\ is a leaf of a connected component of , which consists of 
the path 7^; . 

(IV) (2) If £\^j^i € then 7^1 is a loopless oriented path ending at a vertex of one of 
the cycles of C. The vertex £\ is smaller than the 2ki first vertices of the path, 
but cannot be compared to vertices of the cycle of C. By construction of the 
orientation of MqUM, see Section [2. 11 the orientation of the cycle is compatible 
with that of the edge (11,^2), where ii is the smallest vertex of the cycle and 
12 is its partner in Mq. The vertex l\ is a leaf of a connected component of Fj, 
which is a unicycle with 7^1 as unique branch and a cycle of C as cycle. 

(IV) (3) when £ {£{)', ■ ■ ■ ,4,' (4,)'}' ^^en 7^, contains a loop of length > 3. 
Recall that configurations with odd cycles cancel because of the skew-symmetry 
of the matrix, so that we only consider configurations where iki+i = for some 
s G {1, • • • , fcj}. In this case, the part of the path 7£i to the loop has even length, 
is alternating and start with an edge of Mq. The loop has even length > 4, is 
alternating and its orientation is compatible with the orientation of the edge 
(il The vertex i\ is smaller than all vertices of the path to the cycle and 
smaller than all vertices of the cycle. 

If ^k +i ^1? then i\ is a leaf of a connected component of Fi which is a unicycle 
with a unique branch, consisting of the path 7^* . 

If ^k +i ~ ^1' then i\ is the smallest vertex of a connected component of Fi 
which is a cycle, consisting of the path 7^1 . 

• If 4^_|_i € v"''^'^ ' '^^1 then, the path 7^* attaches itself to a connected component 
of -Fj-i, this can only occur when i > 2, and one of the following happens. 
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(IV) (4) 



The path 7^* attaches itself to a leaf of that is ik^+i = ^1 for some 

t E {!,■ ■ ■ ,i — 1}. Then 7^^ is a loopless oriented path from i\ to i\. The 
vertex £\ is smaller than the 2ki following ones, but greater than than i\. 
Indeed l\ is the starting point of a previous step of the algorithm. This allows 

to identify the ending vertex of the path 7^^ . The vertex i\ is a leaf of Fi . 



(IV) (5) The path 7^1 creates a new branch of the component. Then 7^1 is a loopless 
oriented path. The vertex i\ is smaller than the 2fcj following ones, but we have 
no a priori information on the last vertex of the path. The last vertex of the 



Lemma 9. The oriented edge configuration Fj is an RCRSF compatible with Mq. 

Proof. By definition of the algorithm, the oriented edge configuration Fj contains as many 
edges as MqUM, that is \ V\ edges. By definition, it contains all edges of Mq, that is |V^|/2 
edges, and by the algorithm no doubled edges of T>, that is \V\/2 edges not in Mq. 

By Point (I) the oriented edge configuration Fj has one outgoing edge at every vertex of 
V, which is equivalent to saying that it is an RCRSF such that edges of each component 
are oriented towards its root, and cycles are oriented in one of the two possible directions. 
It thus remains to show that cycles of unicycles are alternating, and have even length 
> 4. By Point (II), the oriented edge configuration Fj has no doubled edge of V, thus 
if Fj has a cycle, it either comes from Point (IV) (2) meaning that it is a cycle of C 
implying that it is even, alternating and has length > 4; or from Point (IV) (3), when it is 
created by the algorithm. Returning to the description of Point (IV) (3) and recalling that 
the contribution of configurations with odd cycles cancel, we know that it has the same 
properties in this case. □ 

For every i G {1, • • • and for every connected component of Fi which is a cycle C 
created by the algorithm {i.e. not a cycle of the initial superimposition) , denote by mc 
the smallest vertex of C. Define Xj to be: 



If Fi has at least one leaf, let yi be the maximum leaf of Fi, else let yi = —00. 

If both Xi and yi are —00, then Fi has no leaves and only contains cycles of the initial 

superimposition Mq U M. This means that the set V is empty, and that F is the initial 
superimposition Mq U M. This has been excluded here, since the full algorithm doesn't 
even start the opening of edges procedure in this case. Thus max{xj,yi} > —00. 

Lemma 10. For every i G {1, • • • ,j}, the initial vertex i\ of Step i is the maximum of Xi 
and yi- 



path 7^i is nevertheless identified as being a fork. The vertex i\ is a leaf of F^. 
Note that the component of l\ might be a unicycle with a unique branch. If 
this is the case, the branch is the path 7^2 and the cycle was created by Case 
(IV) (3) in a previous step of the algorithm. The vertex i\ is thus larger than 
the smallest vertex of the cycle. 




max{mc : C is a cycle-connected component of Fi, but not of C} if {} 7^ 
—00 else. 
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Proof. By Point (IV) above, the vertex i\ is either a leaf of Fi or the smallest vertex of 
a connected component of Fi which is a cycle created by the algorithm, meaning that it 
is not a cycle of C i.e. not a cycle of the initial superimposition Mq U M. Arguing by 
induction, all leaves and smallest vertices of cycle-components of Fi which are not present 
in C, must be initial vertices of steps i of the algorithm for some i € {!,■■■ Moreover 
by construction, the vertex i\ is larger than all previous initial steps of the algorithm, thus 
proving the lemma. □ 

Properties described in Point(IV) also characterize the path once the initial vertex l\ 
is fixed. This can be summarized in the following lemma. 

Lemma 11. Let i\ be the initial vertex of Step i. Then, 

• Suppose that £\ is a leaf of a connected component of Fi . When the connected com- 
ponent is a unicycle rooted on a cycle created by the algorithm, we assume moreover 
that it contains more than one branch. Then, we are in Cases (IV)(1)(2) or (5) and 
the path is characterized as the subpath of the unique path from i\ to the root of 
the connected component, stopping the first time one visits a vertex which: belongs 
to R or to the cycle of the component; is a fork; is smaller than £\. 

• Suppose that £\ is the leaf of a unicycle of Fi rooted on a cycle created by the algorithm 
and containing a unique branch. If i\ is larger than the smallest vertex of the cycle, 
then we are in Case (IV) (5) and the path is the path from i\ to the cycle, stopping 
when the cycle is reached. Else, if £\ is smaller than the smallest vertex of the cycle, 
we are in Case (IV) (3) and the path 'j^i is the path from £\ to the cycle, followed by 
the cycle, with the orientation specified in (IV) (3). 

• // £\ is the smallest vertex of a connected component of Fi which is a cycle created 
by the algorithm, then we are in Case (IV) (3) and the path 7^1 is the cycle, with the 
orientation specified in (IV) (3). 

Remark 12. If the initial superimposition Mq U M consists of cycles only, that is if the 
set T> is empty, then the output of the full algorithm is -F = Mq U M, which consists of 
alternating cycles of even length > 4. The orientation of cycles is specified in Section 12. H 
and cycles are oriented in one of the possible two directions. In this case also, F is an 
RCRSF compatible with Mq. 

3 Proofs and corollaries 

We now prove Theorem [H Corollary [3] and state and prove the line bundle version of the 
result. 

3.1 Proof of Theorem [1] 

In this section, we prove Theorem [TJ Let Mq be a reference perfect matching of G, and let 
F be an RCRSF compatible with Mq, containing kp unicycles. In Lemma [T^ we suppose 
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that F is an output of the full algorithm and identify 2^^ possible perfect matchings M for 
the initial superimposition Mq U AI. Then, we introduce a partial reverse algorithm used 
to define Condition (C) for RCRSFs compatible with Mq. In Proposition [151 we prove that 
an RCRSF compatible with Mq is an output of the full algorithm if and only if it satisfies 
Condition (C), and if this is the case, it is obtained 2^^ times. The remainder of the proof 
consists in showing that contribution of RCRSFs containing unicycles cancel, and that 
only spanning forests remain with the appropriate weight, thus proving Theorem [TJ 

Let F be an RCRSF compatible with Mq, and let kp denote the number of unicycles it 
contains. If kp ^ 0, we let {Ci, • • • , Ckp} be its set of cycles. For every (ei, • • • ,ekp) € 
{0, 1}'^^, define the edge configuration M*-*^^'"' as follows: 



edges of Mq on branches of F 
edges of Mq on the cycle Cj, when ej = 
^ edge of F \ Mq on the cycle Cj , when ej = 1. 



If kp = 0, then the set of cycles of F is {0}, and we set M^^^'' ' '^'^-f) = Mq. 

Lemma 13. For every RSCRSF F compatible with Mq, and every (ei, • • • , e^) S {0, l}*^^, 
the edge configuration M^^'^'"' '^''f^ is a perfect matching of G. 

Proof. If F contains no unicycles, M^^^'"' '''^f^ = Mq and this is immediate. Suppose 
kp 7^ 0. Since Mq is a perfect matching, and since the restriction of Mq and the restriction 
of M^*^!' " '^^'^ to branches of F are the same, all vertices of F \ {V{Ci), • • • , V{Ckp)} are 
incident to exactly one edge of M^^^'"' '^*=f'^. Moreover, by assumption for every j, the 
cycle Cj is alternating, implying that each vertex of V{Cj) is incident to exactly one edge 
of the restriction of Mq and one edge of the restriction of F\Mq to Cj. As a consequence, 
every vertex of V is incident to exactly one edge of M(^1''"'^'»f), proving that it is a perfect 
matching of G. □ 

Lemma 14. Let F he the superimposition of Mq and of an output of the full algorithm. 
Then, the perfect matching M of the initial superimposition Mq U M , must he equal to 
^{£i,-,efc^) {ei,--- ,£kp) e {0,1}'=^. 

Proof. If F is an output of the full algorithm, then by Lemma and Remark [T^ it is 
an RCRSF compatible with Mq, so that V(ei, • • • ,ekp) G {0, 1}*^^, the perfect matching 
^{eir-- lEfej?) jg defined. Suppose that F is an output of the full algorithm with initial 
superimposition Mq U M, where M is not M^^^'"' ''^''f'> for some (ei, • • • ,ekp) € {0, l}^^ . 
Then, Mq U M contains at least one cycle C which is not Ci, • • • ,Ckp. Returning to the 
definition of the algorithm, we know that cycles present in the initial superimposition are 
also present in the output F. This yields a contradiction since F contains exactly the 
cycles Ci,--- ,Ckp. □ 

Partial reverse algorithm 

Input: an RCRSF F compatible with Mq not consisting of cycles only. 
Initialization : Fi = F. 
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Step i, i > 1 



Let i\ be the largest leaf of Fi , and consider the connected component containing i\ . Start 
from i\ along the unique path joining i\ to the root or the cycle of the component, until 
the first time one of the following vertices is reached: 

• the root vertex if the component is a tree, or the cycle if it is a unicycle; 

• a fork; 

• a vertex which is smaller than the leaf i\. 

This yields a loopless path X^i starting from £\, of length > 1. Let Fj+i = Fi\ X^i . If Fj+i 
is empty or contains cycles only, then stop. Else, go to step i + 1. 

End : since edges are removed at every step and since F contains finitely many edges, the 
algorithm ends in finite time A^. 

Definition 3.1. An RCRSF F compatible with Mq is said to satisfy Condition (C) if 
either F consists of cycles only, or if each of the paths A^, • • • , A^-jv, obtained from the 
partial reverse algorithm has even length and starts from an edge of Mq. 

Proposition 15. Let F he an RCRSF compatible with Mq. 

Then, for every (ei, • • • ,5^^) ^ {^i ^}^^ j ^ ^-^ the superimposition of Mq and of an output 
of the full algorithm, with initial superimposition Mq U M^^^'"' ''^''f\ if and only if F 
satisfies condition (C) . The orientation of cycles of F is specified by the proof 



Proof. Let F be an RCRSF compatible with Mq containing kp unicycles, and denote by 
{Ci,--- ,Ckp} its set of cycles. For every (ei,--- jSkp) ^ {0,1}'^^, the edge configura- 
tion M^^^' " '^'^f'-* is well defined, and by Lemma [THl is a perfect matching. We now fix 
(ei,--- ,ekp) e {0,1}''^. Recall that if kp = 0, then the perfect matching M^*^^'"' ''^''f-) is 
simply Mq. 

In the case where kp ^ 0, (ei, • • • ,ekp) = (!,••• , 1), and F consists of cycles only, the 
superimposition Mq U M^^^'"' '^''f^ consists of cycles only, and F = Mq U M^^^'"' '^^f^ is an 
output of the full algorithm. The orientation of the cycles is specified by the choice of 
orientation of Section 12.11 thus proving Proposition [T5j 

Assume that we are not in the above case. Then F is an output of the algorithm with 
initial superimposition Mq U M^^^'"' '''''f^ if and only if there exists a positive integer j and 

a sequence of paths (7^1 , • • • , 7^^ ) € Fj such that F = Mq U M^^^^]... ,7*/^ • 

Lemmas [TU] and [TT] give a characterization of l\ and 7^1 at every step of the algorithm. 
This allows us to define a full reverse algorithm. 

Full reverse algorithm 

Input: an RCRSF F compatible with Mq, (ei,--- T^kp) £ {0,1}'^^ as above, and the 
corresponding perfect matching M*-^^'"' If kp = 0, then m'^^^'"' '^''f^ = Mq. 

Initialization : Fi = F. 
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Step i, i > 1. 

Since Fi is either Fi or is obtained from -Fj-i by removing edges, the set of cycles of Fi is 
included in the set of cycles {Ci, • • • , C^p} of Fi. 

For every connected component of Fi which is a cycle Ca such that = (meaning that 
Ca is not a cycle of the initial superimposition Mq U M*-"^-^' " '^''f'^), let mc^ be the smallest 
vertex of Ca- Define 

{max{mcc • Ca is a cycle-connected component of Fi, and Ea = 0} if {} 7^ 
—00 else. 

If Fi has at least one leaf, let yi be the maximum leaf, else let yi = —00. Note that by 
assumption, we do not have Xi = yi = —00. We let i\ = max{xi,yi}, and 7^1 be the 
oriented path as characterized in Lemma [TTl 

Let Fj+i = Fi\ 7^1 . If the oriented edge configuration Fj+i is empty, or if it consists of 
cycles of the superimposition Mq U M^^^'"' '^''f^ only, then stop; else go to step i + 1. 

End : since edges are removed at every step and since F contains finitely many edges, the 
algorithm ends in finite time j, for some integer j. 

This defines for every RCRSF F compatible with Mq, a sequence of paths 7«i,-- - ,7„j 

1 *i 

such that F is the union of these paths and of cycles of the initial superimposition Mq U 

a consequence, the oriented edge configuration F satisfies Properties (I), 

(II), (IV) (1) — (5). We are thus left with proving that F satisfies Property (III) if and only 

if it satisfies Condition (C) , i. e. we need to show that the paths (7^1 , • • • , 7«j ) all have even 

1 *i 

length, are alternating and start from an edge of Mq if and only if F satisfies Condition 
(C). 

Observe that initial vertices - ■ ■ of the full reverse algorithm consist of initial 
vertices of the partial reverse algorithm, interlaced with smallest vertices of 

components which are cycles. Indeed, the only difference in the partial reverse algorithm 
is that cycles are not removed, but this does not change the characterization of largest 
leaf. 

If £\ is the smallest vertex of a component of Fi which is a cycle, that is if i\ = Xj, then 
7£i is a cycle Ca G {Ci, ■ ■ ■ ,Ckp} such that = 0. Since F is compatible with Mq, the 
cycle has even length and is alternating. The orientation is fixed by the algorithm and j^i 
always satisfies Property (III). 

If i\ is the largest leaf of Fi, that is if £1 = yi, then in all cases except one, which we treat 
below, the path 7^1 is exactly the path A^^/ of the partial reverse algorithm, for some i' < i. 
Condition (C) says that A™' has even length and starts from an edge of Mq. In order to 

show that this is equivalent to satisfying Property (III), we are left with showing that, by 
construction, the path A™/ is always alternating. Suppose that this is not the case, then 
there are at least two edges of the same kind (either in Mq or not in Mq) which follow 
eachother. This implies that there is a vertex v of the path incident to two edges of the 
same kind. Since Mq is a perfect matching, every vertex is incident to exactly one edge of 
Mq, so that we cannot have two edges of Mq following eachother. Thus these two edges 
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do not belong to Mq. Again, since Mq is a perfect matching, the vertex v is also incident 
to an edge of Mq, implying that v is the end of a branch. By construction of the path 
Xji' , the path must stop at v, implying that one of the two edges is not in A™', yielding a 
contradiction. 

We now treat the last case. If i\ is a leaf of a connected component of Fi which is a 
unicycle rooted on a cycle Ca such that = 0, with a unique branch, and such that £\ is 
smaller than the smallest vertex of the cycle. Then the path is the path A^i' followed 
by the cycle with the appropriate orientation. We have to show that 7^* satisfies Property 

(III) if and only if A™' satisfies Condition (C). By Property (III), we know that the part of 

••1 

jfi stopping when the cycle is reached, which is precisely X^i' , has even length and starts 
from an edge of Mq. This is exactly Condition (C), since by the same argument as above, 

the path A=' is alternating. We conclude by observing that since F is compatible with 

'■I 

Mq, the cycle part of 7^1 is alternating, and starts from an edge of Mq by construction of 
the orientation of the cycle. Thus 7£i satisfies Property (III) if and only if F is compatible 
with Mq and satisfies Condition (C). □ 



We denote by G{Mq) the set of RCRSFs compatible with Mq satisfying Condition (C). 
Let F be an RCRSF of Q(Mq), and let kp be its number of unicycles. If kp 7^ 0, then 
for every {ei,--- ,ekp) e {0,1}'=^', denote by p^^^'- '""f^ the version of F obtained from 
the full algorithm with initial superimposition Mq U M^'^'^'"' '^''f\ with the orientation of 
cycles given by Proposition [T5j If kp = 0, then F is obtained exactly once from the full 
algorithm with initial superimposition Mq U Mq. 

Since MqUM^"^^'"' '^'^f^ has exactly X]!^=i cycles, and since MqUMq has none, we have as a 
consequence of the full algorithm, see Equation that the weight WMoi^^^^' " '^'^^^ \-^o) 
is equal to: 

(_l)Effi£» Yl ae iikp^O 

Y\ CLe if kp = 0. 



Recall that T denotes the output of the full algorithm with initial superimposition MqVJM 
for a fixed reference perfect matching Mq and a generic perfect matching M of G. Since 
we now aim at taking the union over all perfects matchings M, we write T as TMoi^)- 

As a consequence of Proposition [151 we have that UmgAI Tmo{M) is equal to: 

( U U f(^^'-'^'=^mo)u( U ^\^o)- (7) 

{F&g(Mo): fcj-^O} (ei,... ,ej^)e{o,l}'=F {Fee(Mo): fcF=0} 

Returning to the definition of the Pfaffian of Equation ([1]) , using Corollary [8] and Equation 
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([7]) in the last line, we deduce that: 

Pf(A) = ^ WMoiM), (Definition of Equation ([I])) 

= Yl HTm„{M)), (by Corollary E]) 

{FGe;(Mo): fcF^O} (£!,...,£. )6{0,l}*f' {FegiMo):kF=0} 

^ . ' ^ ' ' 

Let us show that (/) is equal to zero. As a consequence of Equation ([HD, it is equal to: 

Observing that the term sgn(cr^^^^^(.i,... .e^^)^^^^^) ngg^(n. ■■ .^fej,)^^^^^^ is independent of 
(ffi, • • • , Ekp), we conclude that: 



(/)= Y Mvo(f^^--'-^^-\mo)) n «o E (-1) 



E M%o(f(^--'^^f\Mo)) n «e)(l-l) 

= 0. 

Thus, 

Pf(A)= 5^ t/;Mo(i^\Mo) 

{FGg(l/o):A:F=0} 

= Y sgn(cJMo(nMo)) n 

{Fee(M„):A:j,=0} eeF\Mo 

The set {F G G{Mq) : kp = 0} consists of RCRSFs compatible with Mq containing no 
unicycles, and satisfying Condition (C). Observing that: 

• the set of RCRSFs compatible with Mq, containing no unicycle is exactly the set of 
spanning forests of compatible with Mq of Section [1.11 of the introduction, 

• in the case of spanning forests, the partial reverse algorithm is exactly the trimming 
algorithm of Section II. H 

• Condition (C) of Definition 11.21 and Condition (C) of Definition 13.11 are the same in 
the case of spanning forests, 

• the permutation o"j\,/g(^\A^fo) obtained from the algorithm is the permutation of Defi- 
nition 11.11 
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we deduce that {F € G{Mo) : kp = 0} = T{Mq), and thus conclude the proof of Theorem 

m 

Example. Let us take the reference matching Mq which followed us throughout the 
paper, and consider the three possible perfect matchings Mi, M2, M3 of the graph G given 
in Figure [2J Figure shows the output of the full algorithm with initial superimposition 
Mo U Ml. Since the superimpositions Mq U M2 and Mq U M3 contain doubled edges 
only, the output of the algorithm with these respective initial superimpositions, are the 
configurations themselves. By the Theorem [H configurations M2 and Mi^2,i have opposite 
weights, and configurations M3 and Mi^3^i as well, so that their contributions cancel in the 
Pfaffian. As a consequence, signed weighted half-spanning trees counted by the Pfaffian 
of the matrix A are those of Figure [10] below. 




Figure 10: Black edges of the above configurations are half-spanning trees counted by the 
Pfaffian of the matrix A. 



3.2 Proof of Corollary [3] 

Let us recall the setting: A^ is a skew-symmetric matrix of size {n + r) x {n + r), whose 
column sum is zero, with n even; A is the matrix obtained from A^ by removing the r 
last lines and columns; and G are the graphs naturally constructed from A^ and A in 
the introduction. 

Recall, see Section 12.11 that the signature of the permutation gm assigned to a perfect 
matching M counted by the Pfaffian of A, depends on the ordering of the two elements 
of pairs involved in the perfect matching, but not on the ordering of the pairs themselves. 
Choosing the signature of gm thus amounts to choosing an orientation of edges of the 
perfect matching M. The Pfaffian of A can thus be written as: 

Pf(^) = sgn(aM) n 

MeM eeM 

where the product is over coefficients corresponding to a choice of orientation of edges of 
M, specifying a choice of permutation um- 

Now, it is a known fact that the determinant of a skew-symmetric matrix is equal to the 
square of the Pfaffian: 

det(^) = ( X] ^^'^i^Mo) W X] ^^"^i^M) n ^e)- 
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As in Section [27T1 for every Mq € A^, we choose the permutation aM using the superim- 
position Mq U M. Equation ([1]) thus yields: 

det{A) = Y: sgn(aM,) J] E sgn(cTMo(M))(-l)'^^*'""*'^' ("l)"^*'"""''^' 11 
As a consequence of Theorem [H this can be rewritten as: 

det(A) = Y sgn{aMo) JI E sgn(crMo(F\A'/o)) H "e)' 

AfoeA^ e<=Mo FgJ"{A/o) eeF\A^o 

= E E (sgn(o-A/o) n "e) sgn(o-Mo(nMo)) n 

A/oeA1 FeT{Mo) eGAfo eGF\A'/o 

where o'Mo{f\A4o) defined in Definition 11.11 

We have not yet chosen the permutation assigned to Mq, we do so now. For every 
F G T{Mq), we chose the orientation of Mq to be the orientation of edges induced by the 
spanning forest F: this is precisely o"jvf(j(i?\7v/o)- Combining the product of coefficients 
over oriented edges in Mq and in F \ Mq yields: 

det(^) = Y E sgn(o-Mo(F\Afo))^ JI 

A/oGAI F£T(Mo) eeF 

= E E n-.. 

A^oGM FgT{Mo) eeF 

thus proving Corollary [3l 
Remark 16. 



1. We now give an intrinsic characterization of Upj^^MJ^iMo), not using reference per- 
fect matchings. 

Consider the trimming algorithm of Section 11.11 applied to general spanning forests 
of (not assuming that they are compatible with a reference perfect matching 
Mq). Since the reference perfect matching is not used in the algorithm, everything 
works out in the same way, and the algorithm yields a sequence of paths A^^ , • • • , . 
This yields the following more general form of Definition 11.21 

Definition 3.2. A spanning forest F of is said to satisfy Condition (C) if each of 
the paths A}, • • • , Xi obtained from the trimming algorithm has even length. Let J- 
denote the set of spanning forests of G^ satisfying Condition (C). 

Lemma 17. 

T= U T{Mq). 

MoeM 

Proof. By definition, we have the following immediate inclusion: Uj\/oGA/f -^("^-^o) 

If Mo and Mq are two distinct perfect matchings of G, then T{Mq) D J^{Mq) = 0. 
Indeed suppose there exists a spanning forest F in the intersection. Then, it must be 
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compatible with Mq and Mg, meaning that it contains ah edges of Mq U Mq. Since 
Mq and Mq are distinct, the superimposition MqUMq must contain a cycle, yielding 
a contradiction with the fact that F is a spanning forest. 

Thus it remains to show that given a spanning forest F satisfying Condition (C) there 
exists a perfect matching Mq such that F E J-'(Mo), meaning that F is compatible 
with Mq and satisfies Condition (C) of Definition 11.21 Let -F be a spanning forest 
satisfying Condition (C), and let A^i, • • • , A^jv the sequence of paths obtained from 
the trimming algorithm. For every i E {1, • • • ,N}, let Mo(A^i) consist of half of 
the edges of A^i such that A^i alternates between edges of Mo(A£i) and edges of 
A^» \ Mq{XpJ, starting from an edge of Mo(A^i). Let Mq = U^iMo(A^»). Then F 
is compatible with Mq and satisfies Condition (C) of Definition 11.21 It remains to 
show that Mq is a perfect matching. The edge configuration consists of \V\/2 edges, 
since by construction, it consists of half of the edges of a spanning forest. Moreover, 
since each of the paths A^i, • • • ,X^n has even length, no vertex is incident to two 
edges of Mq, thus proving that Mq is a perfect matching. □ 

As a consequence, Corollary [3] can be rewritten in the more simple form: 
Corollary 18. 

det(^) =^Ylae. 

2. Let H be the set of cycle coverings of the graph G by cycles of even length: a typical 
element ^ G H is of the form ^ = (Ci, • • • , Ck) for some k. Then, since the matrix A 
is skew-symmetric, the determinant of A is equal to: 

det{A)= n (-i)(n«-+n"0 n (-iK«-e 

5=(Ci,-,Cfc)e5 {i:|C,|>4} ggg: {i:\Ci\=2} 

= E n (-2)(n«e) n 

5=(Ci,-,Cfe)gH {i:|C,|>4} {i:|C,;|=2} 

It is also possible to prove Corollary [18] directly, without passing through the Pfaffian, 
by applying the full algorithm to doubled edges of configurations counted by the 
determinant, and by taking into account all edges instead of half of them. 

3.3 Line-bundle matrix-tree theorem for skew-symmetric matrices 

In the whole of this section, we change notations slightly, and we let ^ be a skew-symmetric 
matrix of size nxn, whose column sum is zero, with n even; G = (V, E) denotes the graph 
associated to the matrix A. 

We now state a line-bundle version of the matrix-tree theorem for skew-symmetric matrices 
of Corollary [3l in the spirit of what is done for the Laplacian matrix in |For93j , [Kenllj , 
but first we need a few definitions. 

A C-bundle is a copy of C associated to each vertex v V . The toted space of the 
bundle is the direct sum W = ®v<^v^v A connection ^ on W \s the choice, for each 
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oriented edge of G of linear isomorphism ipij : Cj — > Cj, with the property that 

V'ij- = i^Ji'i that is, we associate to each oriented edge a non-zero complex number 
ipij such that = ipjl ■ We say that ipij is the parallel transport of the connection 
over the edge The monodromy of the connection around an oriented cycle C is the 

complex number to^ = Yl^^gi^e- 

We consider the matrix constructed from the matrix A and the connection ip: 

A cycle-rooted spanning forest of G, also denoted CRSF, is a an oriented edge configu- 
ration spanning vertices of G such that each connected component is a tree rooted on a 
cycle. In all that follows, we assume that cycles have length > 3. Edges of branches of the 
trees are oriented towards the cycle, and the cycle is oriented in one of the two possible 
directions. 

Consider the partial reverse algorithm of Section 13.11 applied to a general CRSF F. Since 
the reference perfect matching plays no role in this algorithm, everything works out in 
the same way, and the algorithm yields a sequence of paths A^i,-- - ,Xin, whose union 
corresponds to branches of F. 

Definition 3.3. A CRSF of G is said to satisfy Condition (C) if all of the paths A^^ , • • • , A^iv 
obtained from the partial reverse algorithm have even length. Let us denote by G the set 
of CRSFs satisfying Condition (C). Then, 

For a generic CRSF F of G, let us denote by (Ci, • • • , Ck) its cycles. 
Corollary 19. 

detiA^)=j2{ n «0( n nM-c^--=>^])- 

F^g {eebranch(F)} {j:|Ci|is odd} 

■ ( n n«e[2-^c:-^c:'0- 

{i:\Ci\is even}ggf^ 

Proof. We expand the determinant of A^ using cycle decompositions, as we have done for 
the determinant of A in Point 2 of Remark [161 Since the matrix A^ is not skew-symmetric, 
we cannot omit odd cycles, and we let H be the set of cycle decompositions of the graph 
G, that is the set of coverings of the graph by disjoint cycles. A typical element of H can 
be written as ^ = {Ci, • • • , Gk}, for some positive integer k. Then, the determinant of the 
matrix A^ is: 

det(^'^)= Yl n (-l)"'""''(n "-^-+n "-^0 n ((-l)aeV'ea-eV'-e). 

C={Ci,-,Cfc)eH {i:|C,|>3} ggct gggr {i:|C,|=2} 

Using the skew-symmetry of the matrix A and the fact that = ip^^ this yields: 

€=(Ci,-,Cfe)eS {i:|C,| isodd} egc^ 

■ n (-i)(n«^[-c:+-^'])- n iu-i). 

{j:|Ci| iseven>4} {i:|C.|=2} g^c^ 
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Note that in a given covering there is always an even number of odd cycles, since otherwise 
there is no covering of the remaining graph by even cycles. We now fix a partial covering 
of the graph by odd cycles, and sum over coverings of the remaining graph by even cycles. 
Since the contribution of the parallel transport to doubled edges cancels out, and since 
the matrix A has columns summing to zero, we then 'open' doubled edges according to 
the full algorithm, using Remark [T6l Everything works out in the same way, with the 
role of R played by odd cycles. In this case though, because of the parallel transport, 
the contribution of RCRSFs do not cancel, but looking at the proof of Theorem [U we 
now precisely what those are. Summing over all partial coverings by odd cycles yields the 
result. □ 

Remark 20. Theorem 1191 can then be specified in the case of bipartite graphs, in which 
case there are no odd cycles, in the case of planar graphs or of graphs embedded on the 
torus etc. 
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